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ABSTRACT 

A 5.2' X 5.2' region toward the young cluster IC 348 has been imaged in the mihimeter 
continuum at 4.0" x 4.9" resolution with the OVRO interferometer to a RMS noise level 
of 0.75 mJy beam^^ at 98 GHz. The data are used to constrain the circumstellar disk 
masses in a cluster environment at an age of ~ 2 Myr. The mosaic encompasses 95 
known members of the IC 348 cluster with a stellar mass distribution that peaks at 0.2- 
0.5 M0. None of the stars are detected in the millimeter continuum at an intensity level 
of 3(7 or greater. The mean observed flux for the ensemble of 95 stars is 0.22 it 0.08 mJy 
Assuming a dust temperature of 20 K, a mass opacity coefficient of Kq = 0.02 cm^ g~^ at 
1300/xm, and a power law index of /3 = 1 for the particle emissivity, these observations 
imply that the 3a upper limit to the disk mass around any individual star is 0.025 Mq, 
and that the average disk mass is 0.002 it 0.001 M©. The absence of disks with masses in 
excess of 0.025 Mq in IC 348 is different at the ~ 3a confidence level from Taurus, where 
~ 14% of the stars in an optically selected sample have such disk masses. Compared 
with the minimum mass needed to form the planets in our solar system (~ 0.01 Mq), the 
lack of massive disks and the low mean disk mass in IC 348 suggest either that planets 
more massive than a few Jupiter masses will form infrequently around 0.2-0.5 M.q stars 
in IC 348, or that the process to form such planets has significantly depleted the disk 
of small dust grains on time scales less than the cluster age of ~ 2 Myr. 

Subject headings: stars:pre-main-sequence — stars:disks — open clusters and associa- 
tions 



1. Introduction 

High resolution imaging at optical, near-infrared, and millimeter-wave wavelengths have con- 
clusively demonstrated that circumstellar disks are ubiquitous around young solar-mass stars (O'Dell 
&: Wong 1996; Padgett et al. 1999; Dutrey et al. 1996). The long-held notion that planets may fre- 
quently form in these disks has been dramatically confirmed with the increasing number of planets 
detected around nearby main-sequence stars (Marcy, Cochran, & Mayor 2000). While a consensus 
has not yet been reached on how planets form in circumstellar disks, current theories propose that 
gas-giant planets form either as gravitational instabilities in an accretion disk (Boss 2001), or from 
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the more gradual build up of planetesimals to form a rocky core followed by accretion of gaseous 
material (Pollack ct al. 1996). In both class of models, the accretion disk lifetime sets an important 
time scale on the planet formation process by constraining the time available for a gravitational 
instability to develop or the rocky core to form before the gas and dust are dissipated. 

Circumstellar accretion disks around young stars are commonly inferred based on the presence 
of infrared emission in excess of the stellar photosphere. In many stars the excess infrared emission 
from l-100;um can be modeled with a geometrically thin, optically thick accretion disk (Lynden- 
Bell & Pringle 1974), although inner disk holes, flaring, and other refinements to this basic model 
are often needed to reproduce the observations in detail (Adams, Shu, & Lada 1988; Calvet et al. 
1991,92; Chiang & Goldreich 1997). Detailed spectral energy distributions from the near- to far- 
infrared have been compiled for only a few nearby star forming regions (e.g. Strom et al. 1989; 
Wilking, Lada, k. Young 1989), but J-, K-, and L-band imaging studies provide an efficient 
means to search for near-infrared excesses in hundreds of stars (often in clusters) of various ages. 
Recent results indicate that at least 50% of solar mass stars at an age of ~ 1 Myr exhibit a near- 
infrared excess, but that the percentage decreases to < 10% for ages 3-10 Myr (Strom et al. 1989; 
Haisch, Lada, & Lada 2001b). These results are consistent with the notion that the lifetime of 
accretion disks around most solar mass stars is ~ 3-10 Myr. 

Strictly speaking, a near-infrared excess is diagnostic of only hot dust with temperatures of 
~ 1000 K. For a solar mass star, these temperatures are found only within ~0.1 AU of the star. 
Most of the dust mass, and presumably the majority of planet formation, will be located at larger 
radii and have temperatures too cool to radiate in the near-infrared. Further, some disks are thought 
to have inner holes such that stars without near-infrared excesses may still contain substantial disk 
masses. (Although lOjUm observations suggest that the transition from an optically thick to an 
optically thin disk occurs on time scales of ~ 0.3 Myr, in which case near-infrared observations may 
in fact provide an indirect probe of the larger scale disk; Skrutskie et al. 1990, see also Haisch, Lada, 
& Lada 2001a). An additional limitation of near-infrared studies is that the emission is optically 
thick and provides little meaningful constraints on the disk mass (see, e.g., Wood et al. 2002). 
Longer wavelength observations are needed to establish the disk lifetime at large orbital radii and 
to provide measures on how disk masses evolve with time. 

Millimeter and submillimeter continuum emission are the best available tracers of dust in the 
cool, outer disk. The most comprehensive continuum surveys to date have been toward young stars 
in Taurus (Beckwith et al. 1990; Osterloh & Beckwith 1995; Motte & Andre 2001), p Oph (Andre k 
Montmerle 1994; Niirnberger et al. 1998; Motte, Andre, & Neri 1998), Lupus (Niirnberger, Chini, 
k Zinnecker 1997), Chamaeleon I (Henning et al. 1993), and Serpens (Testi k Sargent 1998). These 
observations have shown that at an age of ~ 1 Myr, ~ 20-30% of the stars possess a circumstellar disk 
with a mass greater than the minimum mass of the solar nebula (~ 0.01 Mq; Weidenschilling 1977, 
Hayashi 1981), and that the median disk mass is <O.OO4M0. Further insights into disk evolution 
can be obtained by using rich clusters as a tool to empirically measure the mass of disks as a 
function of stellar mass and age in much the same manner near-infrared studies have used clusters 
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to establish the frequency of inner accretion disks (Haisch, Lada, &; Lada 2001b). Continuum 

surveys of chisters are also important since an appreciable fraction of stars in molecular clouds are 
found in rich clusters (Lada et al. 1991; Carpenter 2000), and the high ultraviolet radiation fields 
produced by 0-stars in clusters may influence disk evolution (Johnstone, Hollenbach, & Bally 1998; 
Storzer & Hollenbach 1999; Scally & Clarke 2001). 

The main observational challenge to studying clusters in the millimeter continuum is that high 
angular resolution is required to resolve the stars and distinguish disk emission from the more 
extended molecular cloud. The only rich cluster surveyed to date by interferometric techniques 
has been the Orion Nebula Cluster (Mundy, Looney, k. Lada 1995; Bally et al. 1998). However, 
these observations were 10-50 times less sensitive to dust mass than the observations of Taurus and 
p Oph and did not place stringent limits on the amount of circumstellar mass around individual 
stars. Given the limitations in inferring disk evolution from near-infrared excesses, it is important to 
conduct millimeter continuum observations of additional star forming regions in order to understand 
the evolution of disk masses at young stellar ages when planets are in their formative stages. 

To further measure the circumstellar disk masses in different star forming environments, I 
obtained A3mm continuum observations of the young cluster IC 348 with the OVRO millimeter 
wave interferometer. Based on statistical analysis of near-infrared star counts, the IC 348 cluster 
contains ~400 stars distributed over a ~20' x 20' region (1.9 pc x 1.9 pc at the assumed distance 
of 320 pc; de Zeeuw et al. 1999; Herbig 1998), with about half the stars located within the central 
7' X 7' region (Lada &: Lada 1995; Carpenter 2000). To date, ~200 stars and brown dwarfs have 
been individually identified as likely cluster members based on their spectroscopic and photometric 
properties (Herbig 1998; Luhman et al. 1998; Luhman 1999; Najita, Tiede, & Carr 2000), and 115- 
173 likely members have been detected with x-rays (Preibisch & Zinnecker 2001). A cluster age 
of ~ 2 Myr (Luhman 1999; see also Herbig 1998) has been inferred by placing the low mass stars 
on the HR diagram and using theoretical prc-niain-sequcncc evolutionary tracks from D'Antona 
&; Mazzitelli (1997,98). Near-infrared imaging surveys of the cluster show that ~65% of the stars 
contain a K—L excess characteristic of an optically thick accretion disk (Lada k, Lada 1995; Haisch, 
Lada, k. Lada 2001a). By mosaicking the central 5.2' x 5.2' region (0.48 pc x 0.48 pc) of the IC 348 
cluster in the millimeter continuum, I investigate the circumstellar disk masses in 95 known cluster 
members. 

The OVRO observations of the IC 348 cluster are described in Section 2. The stellar prop- 
erties of the cluster members within the OVRO mosaic boundaries are reviewed in Section 3, and 
constraints on the circumstellar disk masses are derived in Section 4. Section 5 compares these 
results with disk masses in Taurus and the Orion Nebula Cluster, and discuses the implications for 
planet formation. The conclusions are summarized in Section 6. 
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2. OVRO Observations 

A 5.2' X 5.2' region toward the IC 348 cluster was mosaicked in the A3mm continuum using 
the OVRO miUimeter-wave interferometer between September 2001 and January 2002. Continuum 
data were recorded simultaneously in four, 1 GHz wide continuum channels centered at 96.48 GHz, 
97.98 GHz, 100.98 GHz, and 102.48 GHz. The phase and amplitude calibrator was J0336+323 
(a,5 = 03:36:30.10,4-32:18:29.3 J2000), which is located 1.7° from the cluster center. The data 
were calibrated using the OVRO MMA data reduction package (Scovillc ct al. 1993). The adopted 
continuum flux for J0336-I-323, calibrated by observing Neptune and/or Uranus, was 1.60 Jy. The 
RMS dispersion in the measured flux for J0336+323 is 0.10 Jy as computed from observations on 
12 separate nights in which both the calibrator and the planets were observed. 

The mosaic consists of 64 pointing centers arranged in a 8 x 8 grid as shown in Figure 1. The 
pointing centers along a row are separated by 40" (compared to the primary FWHM beam size of 
72"). Adjacent rows are also separated by 40" but with pointing centers shifted by 20" in right 
ascension. The mosaicking sequence consisted of obtaining a 3 minute integration on the phase 
calibrator followed by observing 5 positions in the mosaic for 3 minutes each. This sequence was 
repeated until all 64 positions in the mosaic were observed. The mosaic was observed in the OVRO 
compact, low, equatorial, and high configurations both to increase the UV sampling and to decrease 
the RMS noise. The top few rows in the mosaic were often observed at high airmass with lower 
sensitivity. Additional observations of these rows were obtained to provide more uniform sensitivity 
across the mosaic. 

The 64 individual points were formed into an image using the mosaicking routines in MIRIAD. 

The data were averaged using natural weighting to maximize sensitivity, resulting in a beam size 
of 4.0" X 4.9". Unit gain was achieved over a ~ 5.2' x 5.2' area as shown by the dotted curve in 
Figure 1. The average RMS noise within the unit gain region is 0.75 mJybeam"^ and is spatially 
uniform over the mosaic to within 10%. 

3. The IC 348 Cluster 

Before using the OVRO observations to constrain the circumstellar disk masses in the IC 348 
cluster, I review the stellar properties of the IC 348 cluster members located within the mosaic 
boundaries. The stellar and substellar membership list of the IC 348 chistcr was compiled from 
the spectroscopic observations by Herbig (1998), Luhman ct al. (1998), and Luhman (1999), and 
the narrow band imaging survey by Najita, Ticde, & Carr (2000). To ensure a reliable list of 
cluster members, only sources brighter than Ks=14:.5 were considered cluster members for this 
study, as the field star contamination becomes appreciable at fainter magnitudes (Luhman et al. 
1998; Najita, Tiede, &; Carr 2000). In total, 95 cluster members brighter than iCs=14.5 have been 
identified within the OVRO mosaic. Astrometry for these sources was adopted from 2MASS. The 
astrometric uncertainty (< 0.2") is 20 times less than the synthesized beam size. 
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Stellar masses and ages for 92 of the 95 cluster members have been estimated using the 
D'Antona & Mazzitelli (1997,98) theoretical pre-main-sequence evolutionary tracks and published 
spectroscopic and photometric data as described in Hillenbrand, Meyer, & Carpenter (2002). Three 
sources do not have mass estimates since optical photometry is not available, or in one case, the 
spectral type could not be determined since the source is possibly a Class I object (Luhman et al. 
1998). Figures 2 and 3 show histograms of the inferred stellar masses and ages respectively. (Also 
shown in these figures arc the stellar mass and age histograms for a subsample of the Taurus pop- 
ulation as discussed in Section 5.) The stellar masses range from 0.03 M0 to 4.2 M0 with a broad 
peak in the distribution at about 0.2-0.5 Mq. The mean age of the stellar population is ~2 Myr. 
The properties of this subsample agree well with the properties of the entire IC 348 population 
identified so far (Herbig 1998; Luhman et al. 1998; Luhman 1999). 

To investigate the frequency distribution of inner accretion disks in this subsample of the 
IC 348 cluster, Figure 4 shows the J — H vs. H — Kg and J — H ws. Kg — L color-color diagrams. 
The J, iJ, and Kg photometry are from 2MASS, and the L-band photometry is from Haisch, Lada, 
& Lada (2001a). Only sources with no processing flags in the 2MASS database arc shown. In 
Figure 4, the solid curves represent the locus of main-sequence and giant stars (Bcsscll &: Brett 
1988) and the dashed line is the interstellar reddening vector (Cohen ct al. 1981), where the J — H 
and H — K colors have been transformed into the 2MASS photometric system (Carpenter 2001). 
Of the 82 sources in the J — H ys. H — Kg diagram, 4 have near-infrared colors to the right of 
the reddening vectors that are consistent with an infrared excess. However, for each of these four 
stars, the magnitude of the apparent H — Kg excess is less than the la photometric uncertainties 
and therefore may be attributed to photometric noise. The fraction of stars with an infrared excess 
is higher in the J — H vs. Kg — L diagram since L-band shows larger contrast between the disk 
and photospheric emission. Of the 47 stars with available J, iJ, Kg., and L-band photometry, 18 
(38%) contain an apparent Kg — L excess. Haisch, Lada, & Lada (2001a) do not report L-band 
photometric uncertainties for individual sources. Based on the distribution of stars blueward of the 
un-reddened main-sequence in the J — H vs. Kg — L diagram (a region which cannot be populated 
except from photometric noise and photometric variability between the time of the 2MASS and 
L-band observations) , I estimate that ~ 4 sources may have an apparent Kg — L excess due to 
noise. Therefore, a minimum of ~ 15% of the 95 IC 348 members within the OVRO mosaic have 
near-infrared colors indicative of an optically thick, inner accretion disk. By comparison, Haisch, 
Lada, & Lada (2001a) estimate that the ^65% of the stars in the IC 348 cluster brighter than 
K=12 have a K — L excess. 



4. Constraints on Disk Masses 

Figure 5 shows gray scale and contour maps of the OVRO mosaic. In the contour map, the 
dotted curve shows the unit gain boundary of the mosaic, and the open circles indicate the positions 
of the 95 known members of the IC 348 cluster within the unit gain boundary. The contours begin 
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at 3(7 above the RMS noise of 0.75 mJy beam" with increments of la. Visual inspection of the 

contour map shows that none of the known IC 348 cluster members have been detected in the A3mm 
continuum at the 3a level or greater. Over the entire mosaic, the brightest observed flux is only 
3.9a above the noise, and the frequency distribution of fluxes are consistent with gaussian noise as 
shown by the solid circles in Figure 6. Therefore there are no definitive continuum detections in 
the mosaic. 

Circumstellar disks in Taurus typically have diameters of < 300 AU as measured in the A3mm 
continuum (Dutrey et al. 1996). Assuming that the disks in IC 348 are similar, the disks can 
be considered point sources at the resolution of these observations (1600 AUxl300 AU) when 
computing upper limits to the continuum flux. A histogram of the observed fluxes toward the 95 
cluster members is shown in Figure 6. The dotted line through the histogram shows the expected 
flux distribution for gaussian noise given the observed RMS noise of 0.75 mJy beam~^. The observed 
mean flux toward the cluster members is 0.22 ± 0.08 mJy, where the uncertainty represents the 
standard deviation of the mean of the 95 cluster members. These results indicate that while the 
flux distribution is consistent with gaussian noise, there is a positive bias to the average flux at the 
2.8a confidence level. 

The OVRO continuum observations can be used to place constraints on the circumstellar disk 
masses. In reality, most disks will contain a range of dust temperatures with decreasing tem- 
peratures with increasing radial distance from the star. The data obtained here do not permit 
sophisticated modeling of the disk structure, however, and a single dust temperature was assumed. 
Following Hildebrand (1983), the disk mass (i.e. sum of the dust and gas components) can then be 
estimated as M^isk = ^ f"^^ r, where Ki, = Koiir)^ is the mass opacity coefficient, (3 parame- 

i^u v-t dust ) 

terizes the frequency dependence of the observed flux, T^ust is the dust temperature, and 

By{T(iust) is the Planck function. For consistency with prior studies, I assumed /3=1.0 and = 
0.02 cm^ g^^ at 1300/im (Beckwith et al. 1990). The value of Ko, and consequently the disk masses, 
is uncertain by at least a factor of three (Beckwith et al. 1990). The power law index (5 may be 
as low as ~0.5 on average in circumstellar disks (Beckwith & Sargent 1991; Mannings &; Emerson 
1994). If this lower value of /3 is more appropriate for circumstellar disks, the assumptions adopted 
here will overestimate the disk masses by 50%. 

The appropriate dust temperature for calculating the disk masses was estimated using the 
circumstellar disk models from Beckwith et al. (1990) and Osterloh &; Beckwith (1995). These 
studies fitted radial power-law distributions for the diist temperature and disk mass surface density 
to the observed infrared and millimeter spectral energy distribution for well-studied stars in Taurus- 
Auriga. From the disk mass derived by these model fits, I computed a single component dust 
temperature for each star that reproduces the disk mass given the observed A1.3mm flux. A 
histogram of these 'effective" dust temperatures is shown in Figure 7. (Four stars with effective 
temperatures greater than 100 K are not shown in this figure.) Of the 44 stars in Taurus-Auriga 
with model disk masses from Beckwith et al. (1990) and Osterloh & Beckwith (1995), the median 
dust temperature derived in this manner is ~ 20 K. If this dust temperature is then assumed for 
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all of the Taurus stars, the model-derived dust masses can be reproduced to within a factor of 2 
for 60% of the stars, and within a factor of 3 for 70% of the stars. (For sources with effective dust 
temperatures substantially greater than 20 K, the single component dust model will overestimate 
the disk mass since M^ust cx; T^g^). Assuming that the disk temperatures in the IC 348 cluster 
are not vastly different than in Taurus, adopting a dust temperature of 20 K, while simplistic, 
should reproduce the disk masses within a factor of a few for the majority of the stars. For these 
assumptions, the 3(t upper limit to the circumstellar disk mass around an individual star in the 
IC 348 cluster is 0.025 M0. Treating the 95 cluster members as an ensemble, the observed average 
flux corresponds to an average disk mass of 0.002 it 0.001 Mq. 



5. Discussion 

The limits on the disk masses in the IC 348 cluster can be compared to disk masses in other 
star forming regions as a step toward determining the evolution of disk masses as a function of time 
and environment (e.g. "clustered" vs. "isolated" star formation). The Taurus star forming region 
provides a meaningful comparison since it has been thoroughly studied in the millimeter continuum 
(Beckwith ct al. 1990; Osterloh & Beckwith 1995; see also Skinner, Brown, & Walter 1991 and 
Henning et al. 1998), and the stellar population has an age similar to IC 348 as shown below. 
Since the IC 348 sample has been identified largely based on spectroscopy, only stars in Taurus 
with spectral types and sufficient photometric information to place the stars on the HR diagram 
were used for this comparison. Prom the database of spectral types and photometry compiled from 
the literature by Hillenbrand, Meyer, & Carpenter (2002), 117 stars were identified in Taurus that 
met these criteria and have a published A1.3mm continuum flux measurement. The main source of 
incompleteness in the Taurus sample are deeply embedded IRAS sources which are not detected at 
optical wavelengths (Kenyon Sz Hartmann 1995), and recently identified low mass stars and brown 
dwarfs (Briceno et al. 1998,99; Luhman & Rieke 1998) that have not yet been surveyed in the 
millimeter continuum. 

The open histograms in the top panels of Figures 2 and 3 show the mass and age distributions 
respectively for the Taurus comparison sample. (The bottom panels show the analogous distribu- 
tions for the IC 348 sample as already discussed.) The mean stellar mass of the Taurus sample 
is logio (mass/M©) = —0.33 with a dispersion of 0.27. By comparison, the IC 348 sample has a 
mean logio stellar mass of —0.53 with a dispersion of 0.44. Therefore, the stellar mass distribution 
in IC 348 is skewed toward smaller stellar masses compared to the Taurus sample. Similarly, the 
mean stellar age of the Taurus sample is logio (age/yr) = 6.0 with dispersion of 0.6, while the mean 
logio age and dispersion for IC 348 are 6.3 and 0.6 respectively. Thus the stellar ages for IC 348 are 
skewed toward older values than stars in Taurus. However, the uncertainty in the distances to the 
star forming regions, which is at least ~ 10%, (Kenyon, Dobrzycka, & Hartmann 1994; de Zeeuw 
et al. 1999), will produce a ~0.09 dex uncertainty in the mean age (Hartmann 2001). At best then 
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the difference in the mean ages between Taurus and IC 348 are significant at the ~ 2cr confidence 
level. 

Disk masses for stars in Taurus were re-computed from published A 1.3mm fluxes and using the 
isothermal dust model adopted here for consistency. The hatched histograms in Figures 2 and 3 
show the distribution of stellar masses and ages for stars in Taurus that have disk masses in excess 
of 0.025 M0, which is the 3a upper limit to the disk masses toward individual stars in IC 348. 
Figure 2 shows that in this Taurus sample, such disks masses are found only around stars more 
massive than 0.27 M0. While the lack of such massive disks around lower mass stars may simply be 
due to the relative lack of sources, only stars more massive than 0.27 are used in the following 
comparison between Taurus and IC 348. In the Taurus sample, 14 of the 99 stars with stellar 
masses greater than 0.27 Mq have disk masses in excess of 0.025 Mq. In IC 348, none of the 48 
stars in the same stellar mass range have such disk masses. Using the two-tailed Fisher Exact Test, 
the probability that the frequency distribution of massive disks in IC 348 and Taurus have been 
drawn from the same parent population is 0.0028. This probability is independent of the assumed 
values of Kg and Tdust used to compute the disk masses to the extent that the disk properties in the 
two regions are similar. However, if the vahic of (3 is as low as 0.5 on average (Beckwith & Sargent 
1991; Mannings &: Emerson 1994) as opposed to the assumed value of (3=1, then the disk masses in 
IC 348 will be systematically overestimated compared to Taurus, and the difference in the frequency 
of massive disks will be enhanced. Thus despite the similar stellar ages of the star forming regions, 
the IC 348 cluster lacks the massive disks found in the Taurus at the ~ 3a confidence level. The 
cause of this apparent difference cannot be due to an photoevaporation of the disks in IC 348 from 
a high ultraviolet radiation field (e.g. Johnstone, Hollenbach, &: Bally 1998) since the spectral type 
of the most massive cluster member in IC 348 is B5 (Luhman et al. 1998). Whether or not the 
difference in disk masses can be attributed to the affects of a "cluster" vs. "isolated" environment 
remains speculative given the limited sample of star forming regions, although continuum surveys 
of additional clusters are possible with current interferometers. 

The IC 348 results can also be compared to the disk masses in Orion Nebula Cluster using the 
86 GHz continuum survey by Mundy, Looney, &: Lada (1995). Again using a 20 K dust temperature 
for consistency with the assumptions adopted here, the 3a upper limit to individual disk masses in 
Orion is 0.17 Mq. Mundy, Looney, & Lada (1995) derived a 95% upper limit to the characteristic 
flux at 86 GHz of 1.1 mJy, which corresponds to a disk mass of 0.03 Mq. Thus the results from 
Mundy, Looney, & Lada (1995) exclude the presence of massive disks in Orion, but do not set as 
stringent limits as the IC 348 observations for disk masses in cluster environments. Similarly, Bally 
et al. (1998) derived a 3a upper limit of 0.015 Mq to the disk mass around 5 proplyds in the Orion 
Nebula cluster for an assumed dust temperature of 50 K. If a 20 K dust temperature is instead 
adopted, the corresponding upper limit is 0.047 Mq. 

Finally, the IC 348 results can be contrasted with the minimum disk mass needed to form the 
Solar System. Summing the mass contained in the planets and assuming an interstellar gas to dust 
ratio, the minimum mass of the primitive solar nebula is ~ 0.01 Mq (Weidenschilling 1977; Hayashi 
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1981). This minimum mass reflects primarily the mass needed to form a Jupiter- mass planet. The 
actual mass of the primitive solar nebula may of course be higher depending on the efficiency in 
converting dust and gas in the primitive solar nebula into planets. The upper limit to the typical 
disk mass in IC 348 is comparable then to the minimum mass of the solar nebula, and the average 
mass is ~ 5 times lower than minimum mass of the solar nebula. Given that the disk masses are 
uncertain by at least a factor of 3 due the mass opacity coefficient (kq), these results suggest either 
that the formation of planets with masses greater than a few Jupiter masses is relatively rare around 
0.2-0.5 M0 stars in IC 348, or that the formation of gas-giant planets has significantly depleted the 
disk of small dust grains on a time scale of < 2 Myr, the mean age of the IC 348 cluster. 



6. Summary 

1 present OVRO observations of a 5.2' x 5.2' region toward the IC 348 stellar cluster that has 
been imaged in the 98 GHz continuum to a RMS noise level of 0.75mJybeam~^ at 4.0" x 4.9" 
resolution. A total of 95 known members of the IC 348 cluster brighter than ^^=14.5 are located 
within the mosaicked region. The stellar masses range from 0.03 Mq to 4.2 Mq with a peak in the 
distribution between 0.2 Mq and 0.5 Mq. At least 15% of these stars are surrounded by an optically 
thick, inner accretion disk as evidenced by a Kg — L excess (see Haisch, Lada, & Lada 2001a). The 
OVRO observations arc used to place constraints on the circumstellar disk masses and provide a 
snapshot of disk evolution in a cluster environment at an age of ^^2 Myr (Luhman 1999). 

None of the 95 IC 348 members within the OVRO mosaic were detected at the 3a level or 
greater. The mean flux toward the 95 stars is 0.22 it 0.08 mJy, indicating the observed fluxes are 
consistent with random noise although with a slight bias toward positive fluxes. The millimeter- 
wave continuum fluxes were converted into disk masses (the sum of the gas and dust components) 
assuming a dust temperature of 20 K, a mass opacity coefficient of 0.02 cm^ g~^ at 1300/im, and a 
power law index for the particle emissivity of f3=l. For any individual source, the 3a upper limit 
to the circumstellar disk mass is 0.025 Mq. For the 95 sources as an ensemble, the average disk 
mass is 0.002 ± 0.001 Mq. 

The constraints placed on the disk masses in the IC 348 cluster were compared to the disk 
masses around stars in the Taurus molecular cloud. A subsample of the Taurus population was 
identified that has available optical spectroscopy and photometry as well as published millimeter 
continuum observations. The stellar mass range was restricted to stars with masses in excess of 
0.27 Mq since that is the lowest mass star in the Taurus sample that has an observed disk mass in 
excess of 0.025 Mq. In the Taurus sample, 14 of the 99 stars with stellar masses >0.27 Mq have disk 
masses in excess of 0.025 Mq, while in IC 348, none of the 48 stars in the same stellar mass range 
have such disk masses. The probability that the frequency distribution of massive disks in IC 348 
and Taurus have been drawn from the same parent population is 0.0028, and are thus different at 
the ~ 3a confidence level. These results suggest that the IC 348 stellar population lacks the massive 
disks that characterizes ~ 14% of the Taurus population despite the fact that the two regions have 
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similar ages (Hillenbrand, Meyer, & Carpenter 2002). Whether or not the differences between the 
frequency of massive disks between Taurus and IC 348 is related to "isolated" versus "clustered" 
star forming environments awaits observations of additional star forming regions. 

The implications of the IC 348 observations for planet formation can be placed in context with 
the minimum mass of the primitive solar nebula ('^O.OIM©; Wcidcnschilling 1977, Hayashi 1981) 
needed to form our Solar System. The IC 348 results suggest that if an optically thick accretion 
disk is the necessary seed to form gas-giant planets, then either few stars in IC 348 with stellar 
masses between 0.2 M0 and 0.5 M0 will form planets more massive than a few Jupiter masses, or 
the process to form such massive planets has significantly depleted the disk of small dust grains by 
the age of the IC 348 cluster, or ~ 2 Myr. 
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Owens Valley Radio Observatory, which is supported by the National Science Foundation through 
grant AST-9981546. This publication makes use of data products from the Two Micron All Sky 
Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and 
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Fig. 1. — Schematic of tlie IC 348 mosaic made with OVRO. The sohd squares mark the 64 pointing 
centers used to create the mosaic. The pointing centers are separated by 40" along a row, with 
adjacent rows also separated by 40". For comparison, the circle in the lower right corner shows 
the FWHM beam size of a single OVRO antenna (72") at the observed frequency 98 GHz. The 
dotted curve shows the extent of the OVRO mosaic (~5.2' x 5.2') at the unit gain boundary. The 
synthesized beam size is 4.0" x 4.9". 
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Fig. 2. — Bottom: Histogram of the stellar masses for the IC 348 cluster members located within the 
OVRO mosaic. Three probable IC 348 members within the OVRO mosaic did not have sufficient 
data to infer their stellar properties and are not shown in this figure. Top: Histograms of the 
stellar masses for a comparison sample of stars in the Taurus molecular cloud that have available 
optical spectroscopy, photometry, and submillimctcr continuum observations. The open histogram 
represents all stars in the comparison sample, and the hatched histogram shows the stellar mass 
distribution for sources that have a disk mass greater than the 3cr detection limit of 0.025 
for the IC 348 observations. The stellar masses in both IC 348 and Taurus have been inferred 
by placing the stars in an HR diagram using the database compiled by Hillenbrand, Meyer, & 
Carpenter (2002) and the D'Antona & Mazzitelli (1997,98) pre-main-sequence evolutionary tracks. 
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Fig. 3. — Bottom: Histogram of the stellar ages for the IC 348 cluster members within the OVRO 
mosaic. Top: Histograms of the stellar ages for a comparison sample of stars in the Taurus molecular 
cloud as described in Fig. 2. 



1 71 7 I I I I r— I I I I I I f 




J I I 1 I I 1 1 I 1 I I [_l L 



1 2 

H-K K -L 

s s 



Fig. 4. — Left: J — H vs. H — Kg color-color diagram for 82 IC 348 members in the OVRO mosaic 
that have photometry in the 2MASS database without any error flags. Right: J ~ H vs. Kg — L 
color-color diagram for 47 IC 348 members in the OVRO mosaic with 2MASS photometry and L- 
band photometry from Haisch, Lada, &; Lada (2001a). In each panel, the solid curves represent the 
locus of main-sequence and giant stars Bessell Sz Brett (1988) and the dashed line is the interstellar 
reddening vector (Cohen et al. 1981), where the J — H and H — K colors have been transformed into 
the 2MASS photometric system (Carpenter 2001). Four stars exhibit an apparent H — Kg excess, 
but in each instance, the magnitude of the excess is less than the la photometric uncertainties and 
can be attributed to photometric noise. In the J—H vs. Kg — L diagram, 18 stars show an apparent 
Kg — L excess, of which 14 have an excess larger than the estimated photometric uncertainties. 
Therefore, a minimum of 15% of the 95 IC 348 cluster members within the OVRO mosaic contain 
a detectable Kg — L excess indicative of an optically thick circumstellar disk. 
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Fig. 5. — Left: Grayscale image of the A3mm continuum emission toward IC 348. Darker regions 
represent bright intensities. Right Contour map of the OVRO mosaic. Contours begin at Scr above 
the RMS noise of 0.75 mJybeam"-*^ with increments of la. The dotted boundary that encompasses 

the image shows the unit gain boundary of the mosaic (see Fig. 1). Open circles represent 95 
probable members of the IC 348 cluster within the field of view of the OVRO mosaic that have 
been identified from spectroscopy (Herbig 1998; Luhman et al. 1998; Luhman 1999) and narrow 
band imaging (Najita, Tiede, & Carr 2000). This figure shows that none of the known IC 348 stars 
were detected in the A3mm continuum at the 3cr noise level or greater. 
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Fig. 6. — Frequency distribution of the observed flux densities at 98 GHz over the entire OVRO 
mosaic (sohd circles) and toward the 95 IC 348 members (histogram). The frequency distribution 
indicated by the solid circles have been arbitrarily scaled by 1/600. The dotted curves through 
the circles and through the histogram show the expected distribution of fluxes in the two samples 
for gaussian noise with a dispersion of 0.75mJybcam~^, which is the average noise in the OVRO 
mosaic. This figure shows that the fluxes over the entire mosaic are consistent with gaussian noise. 
The mean flux observed toward the 95 IC 348 members is 0.22 ± 0.08 mJy, where the uncertainty 
is the standard deviation of the mean. 
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Fig. 7. — Histogram of the effective dust temperatures of circumstellar disks in Taurus- Auriga. The 
effective dust temperature represents the dust temperature that, when combined with the observed 
Al.Smm fluxes, reproduces the disk masses derived by fitting power-law distributions for the dust 
temperature and mass surface density to the observed spectral energy distribution of young stars 
in Taurus- Auriga (Bcckwith et al. 1990; Osterloh &: Beckwith 1995). Most stars have effective 
dust temperatures of ~ 20 K, which was used to convert the observed OVRO A3mm fluxes to disk 
masses. Four stars have effective dust temperatures greater than 100 K and are not shown in this 
figure. 



